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The isospin-breaking vector meson decay constants are determined from a QCD sum rule analysis 
of the vector current correlator (0\T (l^fVjf ) \0), using a recently proposed implementation of the 
finite energy sum rule approach. The analysis employs the three-loop version of the OPE, and 
two different families of weight functions. It is shown that the requirement of consistency between 
results obtained using these two different weight families leads to a rather good determination of 
the parameter describing the deviation of the D — 6 condensate term in the OPE from its vacuum 
saturation value, and that the ability to determine this value has non-trivial numerical consequences 
to the analysis. The phenomenological relevance of the results to the extraction of the strange quark 
mass and the determination of the 6 th order chiral low energy constant, Q, is also briefly discussed. 



I. INTRODUCTION 

Because the neutral (a = 3, 8) members of the SU (3) f octet of vector currents, = q~/^\-q (with A a the usual Gell- 
Mann matrices), couple to fermions in the Standard Model, it is possible to use experimental data on the spectral 
functions associated with correlators involving these currents to determine certain quantities of phenomenological 
interest. For example, defining the scalar correlators, H ab (q 2 ), by means of 

i J d*x exp( ?ga; ) (0|T (J» J») |0) = (q„q v - q 2 g^) n ab (q 2 ), (1) 

and the corresponding spectral functions, p ab (q 2 ), as usual, by p ab (q 2 ) = ^lmll ab (q 2 ), one finds that (1) integrating 
the difference p 33 (q 2 ) — p 88 {q 2 ) with the weight function occuring naturally (due to kinematics) in the finite energy 
sum rule (FESR) treatment of hadronic r decays [1] produces a sum rule from which one can, in principle, determine 
the running strange quark mass, m s (p) [2], and (2) integrating the same difference p 33 (q 2 ) — p ss (q 2 ) with weight 
function w(s) — 1/s produces a sum rule from which one can extract one of the 6 th order low-energy constants 
(LEC's), Q, appearing in the 6 th order version of the effective chiral Lagrangian [3]. (See Ref. [4] for a discussion of 
chiral perturbation theory (ChPT) and the method of effective chiral Lagrangians in general, Ref. [5] for the form of 
the 0(q e ) terms in the effective Lagrangian in the most general case, and Ref. [6] for both a discussion of the subset 
of these terms surviving when one restricts one's attention to vacuum correlators, and a definition of Q.) 
Of course, J 3 and J® do not couple separately in the Standard Model, but only in the combination 

jr =4+ (2) 

which gives the light quark (u,d,s) part of the electromagnetic (EM) current. Thus, what is measured in e + e~ — > 
hadrons, is not the desired quantities, p 33 and p ss , separately, but the combination 
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p EM {q2) = p 33 (g2) + -^V) + I/V) • (3) 

In the isospin symmetry limit, p 38 would vanish and, since one could then classify the final hadronic states according 
to their G-parity, it would be straightforward to separate the isovector (33) and isoscalar (88) components of the EM 
spectral function. 

In the presence of isospin breaking, however, this process is no longer so straightforward. The most obvious 
experimental signature of the presence of isospin-breaking in e + e~ — > hadrons is the interference shoulder in the 
e + e~ — ► 7r + 7r~ cross-section in the p-ui region [7]. The e+e~ — > ui — > tt + it~ contribution to p EM is clearly, to leading 
order in isospin-breaking, to be associated with p 38 and hence is usually removed explicitly in analyzing the data. This 
removal is accomplished by (1) fitting the parameters of a model for the total e + e~ — > n + ir~ amplitude, consisting 
of p, ui and possible background contributions, to the experimental data, (2) removing the w contribution once the 
fit has been performed and (3) squaring the modulus of the remaining p contribution and identifying this result with 
the p contribution to p 33 [8]. 

While this procedure does remove one source of isospin breaking contamination from the nominal p 33 so extracted, 
it is easy to see that other such contaminations still remain. Indeed, once one allows isospin-breaking, the physical 
p and u> are admixtures of pure isovector and isoscalar states, the size of the admixture of the "wrong" isospin 
component being governed by the scale of isospin breaking. As a consequence, the intermediate p contribution to 
p 38 , for example, does not vanish. In fact, if one denotes the pure isovector p state by p^ and the pure isoscalar us 
state by one expects p contributions to p 38 from two sources: (1) that due to p^-u;^ mixing (a one particle 
reducible contribution, with coupling of the isovector current to the p^ component and the isoscalar current to 
the u/°) component of the p), and (2) that due to the "direct" (one particle irreducible, IP I) coupling of the p^ 
component to the isoscalar current (such a coupling being unavoidable in any hadronic effective Lagrangian based on 
QCD). Thus, removing the contribution due to the intermediate state u> from the e + e~ — ► cross-section, while 

removing part of the p 38 contribution, does not remove it all. One is then left with, not the desired quantity, p 33 , but 
rather with a combination of p 33 and the residual part of p 38 associated with the intermediate p state (plus possible 
additional such contaminations from elsewhere in the spectrum) . Similar isospin-breaking flavor 38 contributions exist 
for e + e~ — > lo — > 37r, complicating the extraction of the isoscalar spectral function. 

Corrections for such isospin breaking effects, which are unavoidable so long as no process exists in which only one 
of the two neutral flavor currents couples, are thus necessary if one wishes to perform phenomenological analyses of 
the type mentioned above. Such corrections would also be important in performing precision tests of CVC, which 
involve a comparison of p 33 and the charged isovector spectral function p^ + \ measured in hadronic r decays (see, for 
example, Ref. [9]). 

It is easy to see that, to be able to make these corrections (at least in the region below s ~ 2 GeV 2 , where the EM 
spectral function is, experimentally, resonance dominated), it is sufficient to determine the isospin-breaking vector 
meson decay constants. Let us first clarify notation. We define the flavor 3 and 8 vector meson decay constants via 

{Q\j;\V{k))=m v F v e ll {k) (4) 
where V = p, w, <j>, ■ ■ ■, e^(fc) is the vector meson polarization vector, and a = 3, 8. Fp , ' and are non-zero 

(8) (3} f3^ 

in the limit of isospin symmetry; Fp , F^ and F^ zero in the absence of isospin breaking. The experimentally 
determined EM decay constants, F EM , are then given by 

Fy M = F 3 + 1=F 8 . (5) 

Thus, for example, the broad p contribution to p EM , usually taken to be associated purely with p 33 , consists not only 
of a flavor 33 contribution proportional to [F 3 ] , but also of a flavor 38 contribution proportional to -^F 3 F 8 . The oj 

contribution to p EM , similarly, contains both a flavor 88 part proportional to [F 8 ] 2 and a flavor 38 part proportional 
to -^F 3 F 8 . The flavor 38 parts, in both cases, are present only due to isospin breaking, and have to be removed 

from the experimental p and w contributions to p EM in order to obtain the corresponding p contribution to p 33 and 
ui contribution to p 88 . 

It is important to stress at this point that the conventional "few-percent" rule-of-thumb for estimating the size 
of isospin-breaking effects, which might lead one to expect such effects to be numerically negligible, is inapplicable 
in the cases involving p 33 (q 2 ) — p 88 (q 2 ) discussed above. This is true for a number of reasons. First, because the 
difference of spectral functions is itself flavor-breaking, the relative importance of isospin breaking is enhanced by a 
factor ~ 3, characteristic of the inverse of the scale of flavor-breaking. Second, the effect of p-u mixing naturally 
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produces corrections for the p contribution to p 33 and uj contribution to p 88 which are opposite in sign; the effects 
therefore add when the difference is taken. Finally, there is a natural numerical enhancement which makes the size of 
the correction needed to remove the p 3S part of the ui contribution to p EM , and hence isolate p 88 , larger than naively 
expected [10]. The latter two points are discussed in somewhat more detail in Section II below. 

In what follows, we evaluate the isospin-breaking vector meson decay constants by performing a QCD sum rule 
analysis of the isospin-breaking vector current correlator II 38 . The vector meson spectral contributions are, in this case, 
proportional to FyFy, so that a determination of this product, in combination with the experimental determination 
of Fy M , given in terms of F v and F v above, allows a separate determination of F v and F v . The rest of the 
paper is organized as follows. In Section II we discuss qualitative expectations for the pattern of isospin-breaking 
corrections based on the structure of the leading (chiral) order terms in the vector meson effective chiral Lagrangian, 
as well as semi-quantitative expectations for their probable scale which can, using this perspective, be obtained from 
experimental data. In Section III, we discuss briefly the form of QCD sum rules employed (a version of FESR), and 
the advantages of this approach. In Section IV, we discuss the input used for the hadronic and OPE sides of the sum 
rules employed, and present our results. Some advantages of the approach, in particular viz a vis the handling of the 
D = 6 terms in the OPE of the 38 correlator, will also be discussed here. Finally, in Section V we summarize, and 
make some brief comments on the phenomenological significance of our results. 

II. CHIRAL CONSTRAINTS AND THE SCALE OF ISOSPIN-BREAKING CORRECTIONS 

ChPT provides both an underlying conceptual framework and systematic procedure [11] for writing down the most 
general effective Lagrangian relevant to a given set of hadronic states which fully incorporates the symmetries of 
QCD and implements the broken symmetries (such as chiral symmetry) with the same pattern of symmetry breaking 
as occurs in QCD. Although the resulting effective Lagrangian, £ e ff, is non-rcnormalizablc, and hence contains a 
infinite number of low-energy constants (LEC's) (coefficients of those terms allowed by the symmetry arguments that 
go into constructing C e ff), it is possible to formulate the theory in such a way that only a finite number of such terms 
appear to a given order in the chiral, or low-energy expansion. (For so-called "heavy" fields, those whose masses are 
non-vanishing in the chiral limit, this requires a re-formulation in terms of velocity-dependent fields [12,13].) The 
leading order terms in this expansion (in which light quark masses, m q , q = u,d, s, count as 0{q 2 ), with q representing 
some soft external momentum), incorporate the leading constraints associated cither with chiral symmetry, or the 
symmetry pattern of its breaking. 

The heavy field implementation of ChPT given in Ref. [13] for the vector mesons and their interactions with the 
members of the pseudo-Goldstone boson pseudoscalar octet provides some useful information about the pattern of 
isospin-breaking mixing in the vector meson sector. One must bear in mind that, in general, mixing in field theory 
is more complicated than in Hamiltonian quantum mechanics. This is because, in an effective field theory, one in 
general has (and cannot, in fact, avoid) momentum-dependent mixing terms. This means that there are off-diagonal 
elements of the wave- function renormalization matrix and, as a consequence, the relation between the original, unmixed 
fields and the final, renormalized and diagonalizcd fields, involves the product of a symmetric matrix and a rotation 
(and hence is not itself a rotation). (See, for example, the discussion of the effects of tt -t] mixing in Ref. [4]; an 
expanded version is also given as part of the discussion of the treatment of the mixed-isospin axial current correlator 
in Ref. [14].) From the heavy vector meson field formulation, however, one sees that the leading (in chiral order) term 
in C e ff generating isospin-breaking mixing involves one power of the quark mass matrix and no derivatives [13], and 
hence produces no off-diagonal contributions to the wave-function renormalization matrix. The leading order mixing 
effect thus results in a physical p and lo basis which is related by a rotation to the original pure isospin p(°), lo^ 
basis, just as in the quantum mechanical case. The "wrong" isospin admixture in the physical p state is thus 
equal in magnitude, but opposite in sign, to that of the p(°) admixture in the physical to state, at this order in the 
chiral expansion. While higher (chiral) order corrections exist, this pattern should remain approximately valid, even 
at higher order. 

Let us now consider the vector meson decay constants. Since the chiral limit is also SU (3) f symmetric, the leading 
order term generating the vector meson decay constants is necessarily S'C/(3)_F-symmetric. When one now considers 
the effects of flavor- and isospin-symmctry breaking (recalling that both are generated by the quark mass matrix, 
and hence both are produced by the same set of terms in the effective Lagrangian), there are two potential sources 
of flavor and isospin breaking in the vector meson decay constants. The first is that associated with higher order 
terms, involving at least one power of the quark mass matrix, coupling the external photon field to the vector meson 
nonet, the second that induced by the leading quark-mass-dependent term, responsible for mixing, discussed above. 
The leading order mixing effect simply reproduces the standard leading order SU(3)f mixing analysis [13], leading 
to near ideal mixing in the vector meson sector. As is well-known, the combination of ideal mixing and neglect of 
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flavor-breaking in the EM couplings of the unmixed states leads to the prediction that the vector meson EM decay 
constants, measured experimentally in V — > e+e~, [15] should be in the proportions : : F^ =3:1: — \/2, 
where the superscript (0) indicates that the couplings refer to the ideally mixed, but isospin pure vector meson states. 
That this prediction is borne out by experiment represents empirical evidence that, despite the potential SU(3)f- 
breaking photon coupling contributions being of the same formal order as effects induced by mixing, the former are 
numerically suppressed relative to the latter. Since flavor breaking and isospin breaking are generated by the same 
terms in the effective Lagrangian, this implies that isospin breaking in the vector meson decay constants should also 
be dominated by mixing effects. 

If we take this point of view then, up to sub- leading corrections, we find, for the physical p and u> decay constants, 
now including isospin breaking and taking into account the relation Fp ^ ~ 3Fi^ , 

F EM = F (0) _ £F (0) „ F (0) ^ _ 

^ M = ^ 0) +eF(°)^F(°)(l + 3e), (6) 
where e is the leading order mixing angle, defined via 

p = p(°)-e^ ), w = w<V+ep<®. (7) 



We note two relevant features of these results: (1) because of the dominance by mixing, the corrections required to 
convert the pure isovector Fp ^ coupling to the experimental F^ M coupling is opposite in sign to that required to 

convert the pure isoscalar F^ coupling to the F^ M and, (2) because of the pattern of ideal mixing and the numerical 
suppression of the isoscalar current relative to the isovector current in J / f M , the magnitude of the correction is a 
factor of 9 larger in the u> than in the p case. This does not mean that, somehow, isospin breaking has become huge 
in the to case, but rather reflects the expected natural suppression of the w coupling relative to that of the p, which 
is confirmed by experiment. 

In general, but particularly in light of the large numerical enhancement just discussed, it is crucial to find a means 
of estimating the isospin-breaking contributions to the vector meson EM decay constants. In view of the discussion 
above, although there exist additional effects beyond those associated with mixing at leading order, one can get a 
rough idea of the size expected for the isospin-breaking vector meson decay constants by considering experimental 
information on p-to mixing, and ignoring all non-mixing effects. Although crude, this estimate provides an additional 
qualitative constraint for our later sum rule analysis. 

In order to obtain the parameter e describing p-ui mixing at leading order, it is sufficient to determine the off-diagonal 
element, H puJ , of the vector meson self-energy matrix. In the past, values for H puJ in the range ~ —4000 MeV 2 have 
been quoted, based on simplified analyses of e + e~ — ► tt + tt~ data in the interference region which effectively assume 
that the one-particle irreducible o/°)7r + 7r~ vertex is zero, even in the presence of isospin breaking. Since effective 
operators which generate such a coupling exist in the vector meson effective Lagrangian, however, this assumption is 
unphysical (in the sense of being incompatible with QCD). Once one includes contributions to the u> — > 7T7t amplitude 
generated both by p-u mixing and the 1PI vertex (whose strength we will denote by g^l-x), the analysis of the 
experimental data is somewhat more complicated but, in principle, allows a separate determination of both H puJ and 

the isospin-breaking ratio of couplings of the isospin pure states G — g^nn / gfJn [16,17]. 

An important feature of the analysis framework developed in Ref. [16] is that, from it, one understands that the 
smallness of previously quoted errors for n pw is an artifact of the unphysical assumption G = 0, and does not survive 
the more general treatment. It is worth outlining the reason why this is the case since, in doing so, it will become clear 
that it is a difficult task to improve the experimental situation sufficiently to really pin down the mixing contribution. 

The contribution of the physical (i.e., mixed-isospin) lu to the amplitude for e + e~ — > tt + tt~ is obtained experimen- 
tally by determining the timelike pion form factor, F 1T (q 2 ), in the interference region and fitting it to a form 



q 2 — m 2 q 2 



+ background (8) 



where m v are the complex pole positions, m v = rhy — imylV, and the fit parameter, (j), is known as the "Orsay 
phase" . The ui contribution in Eq. (8) is generated by the coupling of the physical u> to n + TT~ which, as discussed 
above, has two sources: 1PR (p(°) -uj( ' mixing), and 1PI (associated with the oj^itit vertex). The physical coupling 
is given, in terms of these contributions, by 
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where, as usual, the superscript (0) indicates couplings of the unmixed isospin pure states. In the (physically plausible) 
approximation that one assumes saturation of the imaginary part of H puJ by tttt intermediate states, one finds 

lmU puJ (m 2 p )^-Gm p r p (10) 

and hence, in the narrow p-u interference region, 

n pw ~ fl pLU - \GrhpTp (11) 

where Yl puJ is now real. The mixing angle e is then given by [16] 

iW™ 2 ) 

e = 



ml - m 2 p 



= -izT-zG, (12) 

where 

_ im p r p n P u,(m 2 ) 

z =— 2' 1 = — — ' V 1 ^) 

m£ - m z p m p T p 

One then finds, upon substitution of Eq. (12) into Eq. (9), that 



9u 



G(l-z) + 



tt P ^{m 2 p ) 

m„r, 



»p± p 



9%i ■ (14) 



In many places in the literature, one finds the approximation m 2 — m 2 ~ im^Fp employed. The deviation of the 
quantity z from 1 represents the error made in employing this approximation. Since Re z ~ 1 and Im z is small 
(~ .2 to .3), one might think it rather safe to set z = 1 in the above analysis (this approximation was, in fact, made 
uniformly in analyses previous to the discussion of Ref. [16]). If it were true that this approximation were reliable, 
then the effect of G in Eq. (14) would cancel exactly [18], and the experimental data would determine the real part 
of n pa) in the interference region with the usually quoted errors (H pLU (m 2 ) = — 3844 ± 271 MeV 2 ; see Ref. [17,19,20] 
and earlier references cited therein). 

Unfortunately, it turns out that the approximation is both misleading and unreliable. The reason is that, although z 
is approximately real and near 1, (1 — z) is dominantly imaginary. Since the denominator of the second term in Eq. (14) 
is also dominantly imaginary, the two terms add nearly constructively. Were the phases of these terms to be actually 
identical, of course, it would be impossible to use experimental information to separate them, no matter how precise 
that data. Fortunately, there is a small phase difference which, at least in principle, means that a determination, 
with sufficient accuracy, of both the magnitude, A, and phase, <f>, of the ui contribution to F n , would allow separate 
determination of G and H puJ - From this, one would be able to reconstruct Tl pul and hence determine e. There is, 
however, a practical difficulty, which limits the accuracy obtainable using present experimental data. This is a result 
of the fact that the phase difference of the two terms is not large, so one is rather close to being in the "identical 
phase" situation, for which a separate extraction of the two terms would be impossible. If one takes the updated 
numerical analyses of Ref. [17], for example, one sees that values of U puJ between —4000 and —8000 MeV 2 are allowed 
by present data (with a central value ~ —6800 MeV 2 ), and that, while the central extraction for G is moderately large 
~ .1, G = is only 2^ a distant. A significant improvement in this situation would require a significant reduction of 
the errors in the determination of the Orsay phase. Unfortunately, the prospects for seeing such an improvement in 
the near future are remote, at present. 

Although, with present experimental accuracy, the errors in the determination of e are disappointingly large, we 
can, nonetheless, as explained above, use the range of values obtained in Refs. [16,17] to help set a rough scale for 
the expected size of those corrections required to go from F^f M to F p , and F^ M to F®. Using the central values for 
the four fits given in Table 1 of Ref. [17], one finds that F p is less than F p M by between 0.3 and 3.8% (the former 
corresponding to fixing G = by hand, the latter to the MOW and A solutions contained in the Table of Ref. [17]) 
and F® greater than F^ M by between 2.6 and 24.6%. We will see that the solutions obtained below via the sum rule 
analysis satisfy these rather loose constraints, and, in fact, provide an alternate determination which is to be favored 
since it both has considerably smaller errors and is free of the uncertainties associated with effectively truncating at 
leading chiral order (which goes into the estimate just discussed). 
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III. QCD SUM RULES AND THE CHOICE OF THE FESR METHOD 



As is well-known, the properties of unitarity and analyticity lead to the existence of various (appropriately sub- 
tracted) dispersion relations for typical hadronic correlators, Tl(q 2 ). By the term "dispersion relations" we mean 
here both those relations based on the Cauchy representation theorem, and those based on Cauchy's theorem (the 
vanishing of the integral of an analytic function over a closed contour entirely in the region of analyticity) . The generic 
term "QCD sum rules" describes those versions of these relations obtained by making kinematic restrictions which 
allow one to take advantage of the asymptotic freedom of QCD. In the first case, one simply writes the representation 
theorem for II (q 2 ) with q 2 large and spacelike; in the second, a choice of the integration contour is made so that at 
least part of the contour lies in the region of large spacelike momenta. In either case, the point of the kinematic 
restriction is to take advantage of the fact that, in the region of large spacelike momenta, computational techniques 
based on the operator product expansion (OPE)/perturbative QCD (pQCD) are applicable. 

The generic form of the sum rules generated using the representation theorem is then, up to possible subtractions, 



where the s t h is the lowest physical threshold in the channel in question and the spectral function, p, is defined as 
usual by p = ilmlL The LHS is to be computed using the OPE/pQCD, while the RHS is given in terms of measured 
spectral data and/or some spectral ansatz (involving a limited number of unknown parameters) for the unmeasured 
part of the spectral function. As first pointed out by Shifman, Vainshtein and Zakharov (SVZ) [21], the utility of 
such relations is greatly improved by Borel transformation. The effect of the Borcl transform is (I) to replace the 
weight 1/ (s — q 2 ) in the hadronic spectral integral on the RHS of Eq. (15) with exp (— q 2 /M 2 ), where M, the Borel 
mass, is a parameter of the transformation, (2) to kill all subtraction terms and (3) to create a factorial suppression 
of the contributions of higher dimensional operators on the OPE side of the equation (c/ {Q 2 ) — > c/(n — l)\M 2n ). 
Ideally, on the hadronic side, one would like to choose M as small as possible, in order to suppress the contributions 
from the large-s part of the integral, for which the spectral function will typically be complicated, and difficult to 
model. In contrast, to improve the convergence of the OPE when truncated at operators of relatively low dimension, 
one wants M as large as possible. To be usable, SVZ sum rules thus require a compromise: one must hope to find 
a "stability window" in M, i.e., a region of M values for which both contributions from the complicated part of the 
spectral distribution (on the hadronic side) and from the highest dimension operators (on the OPE side) are not too 
large. Typically, because of the compromise nature of the choice of stability window, this means that neither the 
contributions from the large-s part of the spectrum, nor that from the highest dimension operator retained on the 
OPE side, are negligible [21-24]. 

For the case of sum rules based on Cauchy's theorem, a common choice of integration contour (and the one we will 
employ in the analysis below) is that shown in Fig. 1. The radius, s , of the circular part of the contour is to be taken 
large enough that the OPE, to the order available, is reliable in the spacelike region of the circle. The resulting sum 
rule is then generically of the form 



where w(q 2 ) is any function of q 2 analytic in the region of integration, and C denotes the circular part of the contour, 
traversed counterclockwise (from above to below the cut). The OPE is then to be used on C (see also the discussion 
to follow), while spectral data and/or a spectral ansatz is to be employed in the "hadronic" integral on the RHS. Such 
sum rules are called, generically, finite energy sum rules (FESR's), and have usually been employed with integer-power 
weights (w(s) — s , with k = 0, 1, 2, • • •) [25,26], though the standard theoretical treatment of hadronic r decays, which 
is of the FESR type, involves a more complicated weight, as determined by kinematics [1]. 

It is useful, for the discussion which follows, to recall at this point the distinction between "local" and "semi-local" 
duality. For a typical hadronic correlator, H(q 2 ), as above, the OPE is expected to be reliable, not only for q 2 large 
and spacelike, but also for those q 2 on any circle of sufficiently large radius in the complex g 2 -plane, apart possibly 
from some region whose size is hadronic in scale about the timelike real axis (where the effects of confinement are 
expected to be important) [27]. The term "local duality", in this context, refers to the postulate (underlying, for 
example, all statements giving the number of subtractions required to make a given dispersion relation convergent in 
QCD) that, once one is in a region q 2 ~ s for which the separation between adjacent resonances is small compared 
to the typical resonance width, the region of validity of the OPE on the circle of radius s extends all the way down 
to the real timelike axis. This is equivalent to saying that the hadronic spectral function can then be calculated 
using the OPE. This is what allows, for example, the hadroproduction ratio, R, in e+e~ scattering to be computed 




(15) 




(16) 
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theoretically in the asymptotic region using perturbative methods. Similarly "semi-local duality" refers to the idea 
that, at somewhat lower ("intermediate") scales, where local duality is no longer valid, nonetheless, averaged over 
some range of (timelike) momenta, the mean values given by using either the actual hadronic spectral function or 
the OPE version thereof should be the same. It is important to understand that, empirically, the condition that 
resonance spacing be much smaller than typical resonance widths is crucial to the validity of local duality. Indeed, 
one can test local duality through the application of various FESR's in the case of the isovector vector channel, for 
which the hadronic spectral function is very accurately measured in hadronic r decays [9]. One finds that, even at 
scales as large as to 2 . ~ 3.2 GeV 2 , and even though the experimental spectral function appears rather featureless in 
this region, nonetheless, local duality is rather poorly satisfied [28]. (Note that in this channel, at these scales, the 
separation of resonances is comparable to their widths.) 

In light of the above discussion, we distinguish three (timelike) kinematic regimes: "low" (the region of narrow, 
well- isolated resonances), "intermediate" (the region of the validity of semi- local duality) and "high" (the region of the 
validity of local duality, for which the hadronic spectral function can be reliably obtained using pQCD/OPE methods). 
The distinction is important because of the ways, described above, in which sum rules are practically implemented. 
Typically, one attempts to use known information, contained in the OPE (i.e., a(Q 2 ) — a s (Q 2 )/Tr, and the values of 
various vacuum condensates), to place constraints on the parameters of the (hopefully, physically- motivated) ansatz 
for the corresponding unknown (or not fully-known) spectral function. This attempt will, of course, be successful only 
if (1) one is able to work at scales for which the OPE, truncated at operators of relatively low dimension, and at a 
given perturbative order for the Wilson coefficients of these operators, is well-converged, and (2) one has a qualitative 
form for the hadronic spectral ansatz which is both physically plausible and does not involve a large number of 
unknown parameters. The latter condition can, realistically, be satisfied only if the scale up to which one needs the 
spectral ansatz is not too far into the intermediate region. Once the higher resonances start to overlap, and various 
multiparticle background processes start to become important, the general form of the spectral function will become 
increasingly difficult to anticipate in advance (until, that is, one is at sufficiently large scales that local duality finally 
becomes valid). 

In the SVZ version of QCD sum rules, as explained above, it is only rarely possible to choose the stability window 
for the analysis in such a way as to avoid contributions from the region where either one does not know the qualitative 
form of the spectral function, or where, even if one does, to implement it fully would involve the use of more free 
parameters than could be determined reliably, given the limited amount of information available in the truncated 
OPE [24]. Conventionally, this problem is dealt with by employing a spectral ansatz in which (1) the low-s region is 
assumed to be dominated by one or two low-lying resonance contributions and (2) the intermediate- and high-s region 
is approximated using the local duality version of the spectral function, which one assumes to start at some "continuum 
threshold", s - It is well-known that this form of "continuum ansatz" represents a rather crude approximation, and 
hence can create significant uncertainties in the analysis if the contributions from the continuum part of the spectral 
integral are large for M values in the stability window. Typically, one attempts to minimize this problem by avoiding 
M values for which this is the case, but the cost of doing so is poorer convergence on the OPE side of the sum rule 
and hence increased uncertainties associated with neglect of higher dimension contributions. 

A similar problem exists for the integer-power weighted version of FESR's. One advantage of the FESR approach 
is that, in contrast to the SVZ method, the choice of scale appearing on the OPE side of the sum rule (s ) is not 
constrained by the requirement of working in a stability window. One is, of course, still forced to choose so not too far 
into the intermediate region, since otherwise the spectral ansatz would be too complicated to allow a reliable analysis. 
The possibility of working at such intermediate scales in the FESR approach can, nonetheless, represent a practical 
advantage over the SVZ approach in certain cases. This is true of those channels for which the stability window of the 
SVZ analysis lies at relatively low M {e.g., M ~ 1 GeV 2 , as found for many applications in the literature). In such 
cases, the larger scale of the FESR analysis leads to an improvement of both the convergence by operator dimension 
and convergence by perturbative order of the Wilson coefficients on the OPE side of the sum rule. Unfortunately, 
this advantage is usually more than offset by an increase in the difficulties associated with the use of the local duality 
approximation in the intermediate region. The reason is obvious: with integer-power weights, the region near the 
timclike real axis where the circular part of the contour joins the cut, does not have the exponential suppression 
present for the "continuum" contributions in the SVZ approach. The problems that result can be quantified in the 
case of the isovector vector channel, where the hadronic spectral function is known experimentally. As shown in 
Ref. [28], the errors in integer-power weighted FESR's, even at scales as high as to 2 , can be very large, despite the 
fact that the OPE at this scale is both dominated by the leading (D = 0) perturbative term, and rather rapidly 
converging. (A more detailed discussion of why it is that one cannot judge the degree of validity of the local duality 
approximation simply by looking at the degree of convergence of the OPE in a given region is given in Ref. [28].) 

These problems, encountered in employing the conventional, integer-power-weighted version of FESR's, are, how- 
ever, not intrinsic to the FESR approach. Indeed, there is at least one example of a non-integer-power-weightcd 
FESR that is very well-satisfied: that giving the hadronic r decay widths in terms of an integral over the circle 
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of radius sq = to 2 of the product of the OPE for the isovector vector current correlator and the weight function 
w T (s) = (l — s/m 2 ) (l + 2s/to 2 ) (where the dominant input parameter in the OPE representation is cj(to 2 ), which 
can be taken as obtained by running the value measured at the Z mass down to the r scale). The reason for the 
success of this sum rule is well-known: the juncture of the cut and circular portions of the FESR contour corresponds 
to the edge of hadronic phase space and hence, because of kinematics, the weight function w T (s) has a (double) zero at 
s = m 2 . This zero suppresses contributions from the portion of the circular part of the contour near the real timelike 
axis for which the OPE representation of the correlator is not reliable (at intermediate scales like m 2 ) [1]. This 
immediately suggests that the appropriate way to implement FESR's in other channels is, not with integer-power 
weights, but rather with weights having a zero at s = s . In Ref. [28] it was shown that, in the isovector vector 
channel, where one can check the procedure explicitly, weight functions of either the form 

-M'-iX 1 ^) • ,17 » 

having a single zero at s — s , or 

-w-O-iyO + O- (18) 

having a double zero at s — so, produce extremely well-satisfied FESR's for a wide range of values of sq and the 
continuous parameter, A. In addition, using only the OPE representation, for a range of A and sq, and fitting the 
parameters of a sum-of-resonances ansatz to this representation, results in a very good reconstruction of the hadronic 
spectral function, including a determination of the p decay constant accurate to within a few % [28] . In order to have 
a compact terminology for use in describing the weight families in Eqs. (17) and (18), we will refer to w s (s) and Wd(s) 
as single-pinch and double-pinch weights, respectively. The freedom to vary A plays a role analogous to that of the 
variation of M within the stability window in an SVZ-style analysis, in that it strengthens the sum rule constraints on 
the parameters of the spectral ansatz by working with a range of different weight profiles. An additional advantage, 
at least if one wishes to determine not just the parameters of the lowest resonance in the channel but also those 
associated with higher resonances, is that the weight function can actually be arranged to be larger in the second 
resonance region than in the first. 

In what follows, in light of its success in the isovector vector channel, we will employ the FESR framework, using 
both the single- and double-pinch weight families, to investigate the isospin-breaking vector current correlator, II 38 , 
defined above. As usual in the FESR framework, we will work at scales as high as possible, compatible with the 
constraint of having a tractable and physically sensible spectral ansatz for s < sq- Since little is known about the 
vector meson resonance spectrum beyond the second excited resonance region, and since including even the second 
excited resonance region would lead to a spectral ansatz with more parameters than are generally tractable for the 
present analysis, we are forced to work at scales no higher than ^2.8 GeV 2 . Since the separation of the first and 
second excited vector meson resonance regions is comparable to the resonance widths (the p' and w' lie at 1419 and 
1452 MeV, the p" and u>" at 1723 MeV and 1649 MeV, respectively [15]), it is clear that, at these scales, we are not 
yet in the region of the validity of local duality. This makes the use of the single- and double-pinch families crucial 
to the reliability of the FESR analysis. In order to maintain as good convergence as possible on the OPE side of 
the two sum rule families, while at the same time allowing enough variation in so to get a good determination of the 
parameters of the spectral ansatz, we also restrict our attention to scales, So, greater than 2 GeV 2 . 

IV. DETAILS OF THE ANALYSIS 

Since the general framework to be employed in the analysis has been outlined in the previous section, it remains 
only to discuss the input required on the hadronic and OPE sides of the various sum rules. 

We begin with the hadronic side. We take, as our ansatz for the hadronic spectral function, a sum of resonance 
contributions. For the scales used in the analysis, the resonances present in the region of the hadronic spectral integral 
are the p,u, <j), p' and u)'. (Although the tails of the p" and w" intrude slightly into the hadronic integration region for 
so near 2.8 GeV 2 , their contributions are strongly suppressed by the zeros in the weight functions. We have checked 
that including an effective, combined p"-u" contribution in the spectral ansatz has negligible effect on the extracted 
p, u and (f> spectral strength parameters.) We thus include contributions, written in terms of Brcit-Wigner resonance 
forms, for all these resonances. Because the separation of the p' and to' is much smaller than either of their widths, and 
also to reduce the number of free parameters in the spectral ansatz, we have combined the latter two contributions. 
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The strong overlap of the two resonances would, in any case, prevent one from being able to sensibly extract separate 
p' and lo' strength parameters by means of any sum rule analysis of n 38 . 
The spectral ansatz then has the form 

p 3 V) = ^= [fpS (q 2 m 2 p ) fj (q 2 - m 2 ) + (q 2 m 2 ) + f^S (q 2 ~ m^,)] , (19) 

where 

5 (q 2 - m 2 ) = ^ (20) 

V ' 7T ( q 2 _ m 2 f + m 2 r 2 

(with T the width of the resonance in question). This expression reduces to 5 {q 2 — to 2 ) in the narrow width approxi- 
mation (NWA). The minus sign in front of f u and the factor of l/4\/3 are conventional; inclusion of the former ensures 
that f u and f p become equal in the limit that the spectral contributions in the p-u region are generated entirely by 
leading order p-Lo mixing. For the combined p'-io' contribution we have taken average values for the effective mass 
and width. f p , f u , and / p ' w ' are free parameters, to be determined from the matching of hadronic and OPE sides 
of the single- and double-pinch sum rules, for a range of so, A values. 

A few comments are in order concerning the form of the ansatz above, and the physical meaning of the parameters 
to be extracted from the analysis which follows. 

The first concerns the need for the inclusion of a <j> contribution. Note that the correlator n 38 is very closely related 
to that, If", obtained by dropping the strange part of the hypercharge current from II 38 (the OPE's for the two 
correlators are, in fact, identical to three-loop order). The latter correlator has been studied in a number of earlier SVZ- 
style analyses [29-32]. In the earliest of these, the NWA was employed for all resonances, and no (f> contribution was 
included in the spectral ansatz [29,30]. As pointed out in Ref. [31], however, the existence of significant cancellations 
between the NWA p and lo contributions (which would be exact in the limit of mixing dominance and equality of p 
and lo masses) means that a <j> contribution, even if significantly smaller than the individual p and lo contributions, 
could nonetheless be important. Performing the sum rule analysis with a <p contribution included shows that this 
is, indeed the case [31]. Including the <f> contribution in the spectral ansatz also cured an unphysical feature of the 
solutions obtained earlier, which did not include it [31]. The analysis of Ref. [31], however, still employed the NWA 
for all resonances. 

The second point concerns the need to incorporate the p width into the analysis. Because, again, of the high degree 
of cancellation between the NWA p and lo contributions, it was pointed out that the precise degree of this cancellation 
might well be sensitive to whether or not the difference between the p and lo widths was retained in the spectral 
ansatz [32]. The subsequent analysis of Ref. [32] showed that this is, indeed, the case: the spectral parameters, /y, 
decrease by factors ~ 6 when one employs the physical widths in place of the NWA. 

The third point concerns the interpretation of the higher resonance strength parameters, f<p and f p >u>'- It is, of 
course, very natural to take the spectral function to be resonance dominated. Moreover, the near-threshold region 
of the spectral function has been computed to two loops in ChPT [33], and one can see from this result that the 
corresponding low-s background contribution to the relevant spectral integrals is tiny compared to that from the 
p-uj region. The case of the 4>, however, is less clear, since background contributions above the p-ui region are not 
amenable to as reliable estimates, and hence might not be similarly negligible. In the ansatz as written, such physical 
contributions, if present, could only be mocked up (approximately) by additional effective contributions to the 4> 
and p 1 -lo' strengths. Thus, one must use some caution in interpreting, for example, the extracted in terms of the 
physical resonance parameters FjJ and F® - some portion of could actually correspond to an averaged version of 
background contributions in the region between p-ui and p'-io'. The quality of the agreement between the hadronic 
and OPE sides of our sum rules is, however, post facto evidence in favor of resonance dominance, and hence also in 
favor of the possibility of interpreting in terms of F 3 and F^. 

Let us turn, then, to the input on the OPE side of the sum rules. We will discuss the contributions, in turn, by 
operator dimension. 

Since the correlator in question is isospin-breaking, the only dimension D = contribution to II 38 is electromagnetic 
(we adhere, here, to common usage, according to which the leading mass-dependent perturbative terms are labelled 
D = 2). We retain only the leading order (2-loop) graph in this case. 

The D = 2 contributions are dominated by the strong interaction terms proportional to (m<j — m u ) 2 . To 3- loop 
order, the results for these terms follow from the 3-loop expressions for the correlator involving a flavor-non-diagonal 
current and its conjugate [34], since the perturbative contributions involving two quark loops and a purely gluonic 
intermediate state (present for flavor diagonal currents but not for flavor-non-diagonal currents) do not enter until 
4-loop order. The resulting expressions are given in the Appendix. To evaluate them, we require the running masses, 
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m(Q 2 ), and running strong coupling, a s (Q 2 ). These can be obtained once the values are determined at any fixed 
scale, jLto- Since the 4- loop 7 [35] and (3 [36] functions for QCD are now known, we have employed these when running 
the masses and coupling (explicitly, we solve the RG equations exactly, using the truncated 4-loop 7 and (3 functions 
as input). 

As input for the running coupling, we take [i = m T and use the latest (1998) value for a s (m 2 ) obtained by the 
ALEPH Collaboration in their analysis of non-strange hadronic r decays [37] . (The analysis of the strange decays 
employed previous theoretical results for the D = 2 terms, proportional to (m s — m u ) 2 , which turn out to be in error 
[38-40]; the value obtained in the global ALEPH analysis must, therefore, be excluded.) 

The situation with the light quark mass difference 5m = (m d — m u ) is somewhat more complicated. We first write 

5m = [ md W " ) (m d + m u ) = r(m d + m u ) . (21) 
\m d + m u ) 

The isospin-breaking mass ratio, r, is known, from a number of ChPT analyses, to be r — 0.288 ± 0.037 [41], which 
would allow one to determine 5m if m d + m u were known. The most recent determination of m d + m u is that based on 
an integer-power-weighted FESR analysis of the isovector pseudoscalar channel [26,42]. In this analysis, the pion pole 
contribution to the spectral function is known experimentally but the continuum contribution is not. The authors 
of Refs. [26,42], therefore, constructed an ansatz for the unmeasured continuum contribution. It turns out that the 
continuum portion of the resulting model spectral function provides roughly 3/4 of the contribution to the extracted 
value of (m d + m u ) 2 . Unfortunately, it has recently been shown, using the FESR framework discussed above, that 
this continuum ansatz is unphysical [28], so one cannot employ the values of Refs. [26,42]. 

If m s were known (at some scale) , then one could straightforwardly determine m d + m u (at that same scale) using 
the known (scale-independent) ratio of masses, r s — 2m s /(m d + m u ) = 24.4 ± 1.5, obtained by Leutwyler [41] using 
ChPT. Unfortunately, the situation is also somewhat complicated for m s . A number of recent analyses produce 
values of m s (l GcV 2 ) (in the MS scheme) ranging from <~ 110 McV to ~ 210 MeV, often with rather large errors 
[43-47,38,39,48]. Because the analyses based either on flavor breaking in hadronic r decays [38,39] or Narison's r- 
decay-like sum rule for II 33 — II 88 [46,47] have rather large errors resulting from experimental uncertainties which 
are unlikely to be significantly improved in the near future, the most favorable approach would appear to be that 
based on various sum rule treatments of the strange scalar channel, where the dominant Ktt part of the spectral 
function is in principle determined, via the Omnes representation of the timelike scalar Ktt form factor, in terms 
of experimental Ktt phase shifts and K e ^ data [43]. The most recent analyses of this channel [45,48] employ the 
SVZ framework, and produce values m s (l GeV 2 ) = 125 - 160 McV [45], and 160 ± 30 MeV [48]. (The same low-s 
part of the spectral function is used in both analyses; the only difference between the two lies in the treatment 
of the "continuum". The results of Ref. [48], in addition, show no stability window for m s ). Preliminary work 
using the FESR framework discussed above indicates that residual errors associated with the use of the local duality 
approximation in the continuum region remain, for this channel, when one uses the SVZ approach. (See, e.g., the 
results of Ref. [28]. From these one can see (1) that using the central values for the parameters describing the fit to 
the Ktt phases from Refs. [43,45], together with the central value from the m s range from Ref. [45], one obtains rather 
poorly satisfied families of FESR's, and (2) that using the spectral function of Refs. [43,45,48], again with central 
values for the fit parameters, the FESR analysis, in fact, produces m s values larger by ~ 20 MeV than those obtained 
in the analysis of Ref. [45]). Although work on the extraction of m s is still in progress [49], we conclude already from 
the preliminary results noted above, that m s (l GeV 2 ) ~ 165 MeV, probably with errors <~ ±20 McV or less. The 
ChPT ratio then produces (m d + m u )(l GcV ) ~ 13.5 MeV, with errors ~ ±2 McV. For any value in this range it 
turns out that the D — 2 contributions are at the <~ 15% (or less) level of the D = 4 contributions, and the resulting 
errors lead, therefore, to very small (% level) uncertainties in the final results. Since these uncertainties are much 
smaller than those generated by the uncertainty in the isospin-breaking mass ratio r, we have employed the central 
value (m d + m u )(l GeV 2 ) = 13.5 MeV, and retained only the uncertainty in r, in the analysis which follows. 

The D — 4 contributions are much more straightforward. Although in principle both those D = 4 terms proportional 
to the isospin-breaking mass difference, 5m, and those proportional to the isospin-breaking condensate difference, 
< dd > — < uu >, appear in the OPE of II 38 , the latter are numerically tiny compared to the former. The dominant 
D = 4 contribution can then be written in terms of r and the combination (m d + m u ) < qq >, which we can take 
from the GMO relation 

(m d + m u ) < qq >= -mlf 2 . (22) 

The dominant uncertainties for the D = 4 terms thus result from those on r. 

The phenomenological situation is not so favorable in the case of the D = 6 condensates. Usually, in the absence 
of pre-existing determinations of the relevant condensates, one makes estimates based on the vacuum saturation 



10 



approximation (VSA). It is well known that, in situations where it has been possible to perform phenomenological 
checks by extracting the total D = 6 contribution from data, the VSA has proven to significantly underestimate these 
contributions [50]. Usually one simply replaces the factor a s < qq > 2 , which is produced by the VSA, by an effective 
scale-independent factor, written p'a s < qq > 2 . The parameter, p', then represents the deviation from the VSA. 
Ideally, it should either be possible to determine p' from data, or the D = 6 contributions should be small, for the 
sum rule in question. In our case, neither of these conditions holds. In particular, because we are forced to work at 
scales as low as 2 GeV 2 in order to constrain the spectral parameters, the D — 6 contributions can, for so ~ 2 GeV 2 , 
and certain values of A employed in our analysis, approach ~ 40% of the leading D = 4 term. Fortunately, it turns 
out, as we will see explicitly below, that by working with both the single- and double-pinch weight families, we can 
actually obtain a rather good determination of the D — 6 contribution to the correlator (albeit it as a function of r) 
by insisting on the consistency of the results obtained from the two different sum rule families. 

In the Appendix, it is shown that the VSA leads to an expression for the D — 6 contribution to II 38 proportional 

to 



a s ((dd) 2 - (uu) 2 ) = j(a s (qq) 2 ) , (23) 
where (qq) is the average of the u and d condensates, and 

(uu) 

describes isospin-brcaking in the light quark condensates. In order to compare the deviation from the VSA in the 
isospin-breaking channel with that in the analogous isospin-conserving isovector vector (ab = 33) channel, we write 
the re-scaled version of the RHS of Eq.(23) in the form 

p'a s ((dd) 2 - (uu) 2 ) = p red 7 (pa s (qq) 2 ) , (25) 

where p is the parameter describing the deviation from the VSA in the 33 channel, and one has, phenomenologically, 
[51] 

pa s (qq) 2 = (5.8 ± 0.9) x 10~ 4 GeV 6 . (26) 

With this definition, p re d reduces to 1 in the limit that the deviation from the VSA is the same in the 33 and 38 
channels. 

The consistency procedure for fixing the D — 6 contribution to II 38 , together with the phenomenological input 
of Eq. (26), of course, determines only the product p re dl- In presenting our results for p re d below, we have taken 
7 ~ —0.008, which represents an average of the previous determinations listed in Ref. [30], bar one. (We omit the value 
based on an analysis of baryon splittings because it implies (via the 1-loop ChPT relations between flavor-breaking 
and isospin-breaking in the light quark condensate [4]) (ss) / (uu) > 1, which appears unphysical). We will discuss the 
determination of p re d in more detail below when we present the results of the analysis. 

The last point in need of discussion concerns the way in which we handle the integrals on the OPE side of the 
various FESR's. Two options exist in the literature. The first, sometimes called the "fixed order expansion", involves 
first expanding a s (Q 2 ) and the mass factors, generically m(Q 2 ), in terms of a s (so) and m(so). The coefficients of the 
perturbative expansions in powers of a s (so) are then polynomials in £n(s/s ) [52], and the desired contour integrals can 
thus be written in terms of elementary integrals involving logarithms and powers of s. The integrated OPE expressions 
which result involve m(so) 2 , multiplied by a power series in a s (so)- There is, of course, in this expression, the usual 
residual dependence on the choice of scale sq for the expansions discussed above, which results from truncating the 
full perturbative series at fixed order. The second alternative, often referred to as "contour improvement", involves 

numerically integrating the factors [ra(Q 2 )] [a s (<5 2 )]"' s l around the circular contour in the s = —Q 2 plane [53]. It 
is known that this has the effect of simultaneously improving the convergence of the perturbative series and reducing 
the residual scale dependence [53,26,42]. As a result, we have evaluated all the integrals on the OPE sides of our sum 
rules using this approach. 

Let us now turn to the results, which are presented in the Table. As explained above, the dominant uncertainty 
is due to that in the ChPT determination of r. We have, therefore, tabulated the results for the range of values 
corresponding to the errors on r quoted by Leutwyler [41]. All results are based on matching the hadronic and OPE 
sides of the two sum rule families for so in the range 2.0 to 2.8 GeV 2 , and with A in the range 2 to 5 for the single-pinch 
case and 3 to 6 in the double-pinch case. The choice of range of A in each case has been made so as to keep the 
convergence of the perturbative series for the D = 2 term under control. It is worth mentioning that the quality of the 
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match between the OPE and hadronic sides which results after the fitting of the spectral parameters is significantly 
better for Leutwyler's central value of r. 

The value of the D = 6 VSA-violating parameter, p rec i, given in the Table, is determined by requiring that the 
valucds of f p obtained using the single- and double-pinch weight families are the same. The sensitivity of f p to 
variations in p rec i (true also of the other fv), as well as the difference in the p rec i dependence of f p for the single- and 
double-pinch analyses, is shown in Figure 2. The fact that, once p re d has been determined by the requirement of the 
consistency of the two output f p values, all the rest of the spectral parameters, determined using either the singlc- 
or double-pinch weights, also become consistent is strong evidence in favor of the reliability of the analysis. Note 
that (1) the possibility of determining the correction to the VSA for the D = 6 operators, and (2) the inclusion of 
both the D = 2 terms and the 0(a s ,a 2 s ) contributions to the Wilson coefficient of the D = 4 term, are features not 
present in previous analyses of the analogous isospin-breaking IF W correlator. Although the value of p re d, determined 
as just described, depends somewhat on r, this dependence is not strong, and we obtain p rex [ = 1.15 ± 0.15 ± 0.2. 
The first error corresponds to that in Eq. (26), the second to that on r. We see that the violation of the VSA is very 
similar in both the 33 and 38 channels. The importance, in reducing the errors on the determinations of the spectral 
parameters, fv, of being able to determine p re d is also evident from Fig. 2. 

Having determined the D = 6 contributions by self-consistency, the errors on the extracted values of fv are 
determined solely by those on r, and arc <~ 10 — 15%, completely correlated with r. 

Having extracted the parameters fv, it is straightforward to determine the isospin violating decay constants. One 
finds 

F 8 = 2.4 ± 0.3 MeV 
F 3 = -3.4 ±0.4 MeV 

UJ 

Fj = 0.33 ± 0.02 MeV , (27) 
where the errors reflect those on the input isospin-breaking mass ratio, r. 



V. SUMMARY AND DISCUSSION OF PHENOMENOLOGICAL CONSEQUENCES 

A number of useful general observations follow from the analysis above. First, we have found that the violation 
of the VSA for the D = 6 condensates is very similar in the isospin-breaking (38) and isospin-conserving (33) vector 
current channels. Second, we have seen that, although the analysis can be performed successfully for any r in the 
range given by Leutwyler, the central value of that range is preferred, in the sense of giving the best match between 
OPE and hadronic sides of both the single- and double-pinch families of sum rules. Finally, we have demonstrated that 
the FESR method, particularly when implemented using both the single- and double-pinch weight families, is very 
effective, allowing a determination of the p 38 spectral parameters, fv, with rather small errors. These errors (between 
10 and 15% for f p , f u and /^) are a factor of 3 smaller than those obtained in the earlier analysis of Ref. [10] based 
on results of an SVZ analysis of U pul [32]. In order to obtain this level of reduction, the ability to self-consistently 
determine p rec i was crucial. 

Let us now turn to the phenomenological consequences of our results. First note that the corrections required to 
convert the measured contribution of the vector meson, V, to the EM spectral function, p EM , into the corresponding 
contribution to either p 33 (for V = p) or p 88 (for V = w, 4>) are given by the ratios 



F 3 

p 

T?EM 
P 



F 8 



Fl 



^ F EM 



= 0.982 ± 0.0021 



n 2 



= 1.154± 0.017 



= 1.009 ±0.001 , 



(28) 



where the numerical values follow from those in Eq. (27). The size of the deviations of the p and u> corrections from 
1 are reduced by ~ 15 — 20% from those obtained in the earlier analysis [10]; that for the (f> is increased, but remains 
small. In all cases the errors have been reduced by a factor of 3 or more. Note that the first of these corrections is 
the one relevant to precision tests of CVC. Note also that, as claimed above, the corrections given in Eqs. (28), for 
both the p and oj, lie in the corresponding ranges produced by the estimate of Section II. 
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With the results given in Eq. (28), it is now possible to correct the EM data used as input to the inverse moment 
chiral sum rule for the 6 th order LEC, Q. The sum rule is given by [6,54] 



+ 



2tt 2 F 2 



2 / ""K 



(29) 



where \i is the renormalization scale of the effective chiral theory and L r k are the usual renormalized 4 th order LEC's 
of Gasser and Leutwyler [4]. This sum rule was evaluated in Ref. [54] using as input EM data for the isoscalar 
spectral function and both EM and r decay data for the isovector spectral function. The corrections above, required 
for the EM data, were not considered in this analysis. It is not clear, from our reading of the discussion of Ref. [54], 
exactly what the relative weightings of r decay and EM data in the determination of the p° contribution to the LHS 
above actually were. Since, however, r decay data is considerably more precise than electroproduction data, we have 
assumed in what follows that the determination is dominated by r decay data. To the extent that this is true, we 
need only make corrections to the (nominally) isoscalar u and 4> contributions. The result of this exercise is a shift of 
Q(m 2 ) from (3.7 ± 2.0) x 10" 5 to 

Q(m 2 p ) = (2.4 ± 2.0) x 10~ 5 . (30) 

(For reference, making, instead, the somewhat perverse assumption that the determination of the p° contribution was 
dominated by EM data, one would find Q{m 2 ) — (2.0 ± 2.0) x 10~ 5 . The full correction is dominated by that to the 
u contribution. The reason this correction is so much larger than the others has been discussed above.) 

One should bear in mind, in interpreting these results, that (1) there are, in principle, additional corrections to be 
made to the nominal isovector and isoscalar contributions at higher s, and (2) the KK2tt contributions were taken 
to be purely isoscalar in the analysis of Ref. [54] . Because the separations within the higher isovector and isoscalar 
resonance pairs, p'-uj' and p"-u>", are much smaller than the resonance widths, it is not possible to use sum rule 
methods to extract the individual isospin-breaking decay constants of these resonances. As such, we are unable to 
estimate the size of the former corrections. Were the KK2ir states to have a significant isovector component, the 
effect would be to raise Q(m 2 ). 

An alternate method for determining Q is based on the observation that Q occurs not only in the inverse chiral 
moment sum rule above, but also in the 2-loop ChPT expression for n 38 (0) [33]. It is, thus, possible to make an 
independent estimate by using the fitted spectral ansatz for II 38 to compute II 38 (0), assuming negligible contribution 
from the portion of the spectrum above 2.8 GeV 2 . One obtains, from this exercise, 

Q(m 2 ) = (3.3 ± 0.4) x 10~ 5 . (31) 

Since there is no positivity constraint on p 38 (s) one does not know in which direction this result would be changed by 
corrections due to the small higher-s part of the spectral integral. The results of a study of the effect of including two 
combined spectral contributions, one for the p'-uo' and one for the p"-uo" ', however, shows negligible change in Q(m 2 ), 
suggesting that such corrections are unlikely to be numerically significant. Since the two independent determinations 
of Q are completely consistent, within errors, the conclusions that Q{m 2 ) ~ 3 x 10~ 5 is considerably strengthened. 

The last phenomenological application of our results concerns the effect on Narison's r-decay-like sum rule for m s . 
Since a detailed discussion of the way in which one implements the isospin-breaking corrections is given in Ref. [47], 
we report here only the results of employing the improved determinations of the correction factors determined above. 
In doing so we will also take the opportunity to update the input parameters to the analysis of Ref. [47] , employing 
the newer (1998) ALEPH value of a s (m 2 ) [37]. One finds, for example, using r decay data for the isovector input, 
that the average over the values of m s (l GeV 2 ) extracted using the set of scales so = 1-4, 1.5 and 1.6 GeV 2 in the sum 
rule analysis is shifted from 138 MeV to 146 McV (147 MeV if one retains the 1997 ALEPH value of a s as input). 
Unfortunately, the errors on this value associated with uncertainties in the experimental input are still very large, 
<~ ±50 MeV at least, and this uncertainty cannot be appreciably reduced without a significant improvement in the 
accuracy of the determination of the experimental u — > e+e~ and <j> — » e+e~ widths. As such, although the central 
value is brought into better agreement with that discussed above, little more can learned from the Narison sum rule, 
at the present time. 
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APPENDIX: THE OPE FOR IT 8 

The explicit form of the OPE for n 38 , keeping terms only up to dimension six, and to 0(a 2 , m 2 ), can be obtained, 
as explained in the text, from the relevant expressions for the flavor-non-diagonal case given in the literature (see 
Ref. [34] and the paper by Braaten, Narison and Pich (BNP) in Ref. [1]). We list the results by operator dimension. 

Dimension 0: 

The only isospin-breaking contribution at dimension is that due to EM, and is given by [29] 

a 1 



[ n i7£]i) = o 



16tt 3 4V3 



ln(Q ) 



(Al) 



where a is the usual EM coupling. 
Dimension 2: 

The D — 2 term consists of the leading mass-dependent part of the perturbative contribution to the OPE, and 
follows from the expression given in Ref. [34] . One finds 



q 2 [n 38 (Q 2 )] D= 



1 



2tt 2 4V3 



+ 



[(m 2 -m 2 u )(Q 2 



i + 3«(Q 2 ) 



17981 
432 



1045 

~1a 



C(5) a 2 (Q 2 ) 



(A2) 



where a(Q 2 ) = a s (Q 2 )/ir, and ((n) is the Riemann zeta function. Further details on how the running of the coupling 
and the masses is handled can be found in Section IV. 



Dimension 4: 

Our expression for the D — 4 contribution also follows from that given in Ref. [34]. Only the m 4 . and quark 
condensate terms survive once one takes the relevant isospin-breaking difference. The former are numerically tiny 
compared to the latter, and hence have not been written down explicitly. We then find 



4 [n 38 (Q 2 )] D=4 = 



2 ((m u uu) — (rriddd}) 
4\/3 



l + ^(Q 2 ) + y« 2 (Q 2 ) 



(A3) 



The scale- invariant < raqq > difference can be written in terms of the — m u , < dd — uu >, and the averages of 
the u and d quark masses and condensates. Since isospin-breaking in the condensates is much smaller than in the 
masses, the term proportional to md — m u dominates numerically. It can be recast in terms of the isospin-breaking 
quark mass ratio, r, and f n , m^, as explained in the text. 



Dimension 6: 

The 4-quark operators are the dominant operators at dimension 6. Their contribution to II 38 can be obtained 
from the expressions given in the Appendix of BNP [1] . Since lack of phenomcnological information on the various 
condensates forces one to work with the re-scaled version of the VSA, one must, for consistency, drop the terms of 
0(a 2 ) contained there. (See the discussion of this point contained in BNP.) One then finds 

4V3Q 6 [n 38 (Q 2 )] D=6 = -8ir 2 a(Q 2 ) ((u^ l5 T a uu^ l5 T a u) - (d^ l5 T a dd^ l5 T a d)) (A4) 

-^a(Q 2 ) E {(u^T a uq k ^T a q k ) - {d^T a dq kl »T a q k )) , (A5) 
y k 

where T a is an SU (3) generator. Implementing the re-scaled VSA, this expression reduces to 

4V3Q 6 [II 38 (Q 2 )] D=6 = ^L Predl { P a s {qq) 2 ) , (A6) 
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where 7, p and p re d are as denned in the text. 



TABLE I. Results for the D — 6 VSA-violation parameter, pred, and the spectral strength parameters, f p , f w , f$ and f p i w i 
as a function of the isospin-breaking mass ratio, r. The first line, for each value of r, corresponds to the results obtained using 
the single-pinch weight family, the second line to those obtained using the double-pinch family. 



r 


Pred 


/p(x10 3 ) 


/4xio 3 ) 


WxlO 3 ) 


W(xl0 a ) 


0.251 


1.02 


2.3 


1.7 


-0.28 


-0.020 






2.3 


1.7 


-0.28 


-0.020 


0.288 


1.15 


2.6 


2.0 


-0.32 


-0.026 






2.6 


2.0 


-0.32 


-0.026 


0.325 


1.28 


2.9 


2.2 


-0.36 


-0.032 






2.9 


2.2 


-0.36 


-0.032 




FIG. 1. The FESR "Pac-man" contour 
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FIG. 2. The variation of f p with p red for the single- and double-pinch weight families. Results are displayed here for the 
central value r — 0.288. The solid line corresponds to the single-pinch weight analysis, the dashed line to the double-pinch 
analysis. The intersection point determines the value of p re d quoted in the table. 
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